Draft version March 8. 2013 

Preprint typeset using I^T^X style cmulatcapj v. 8/13/10 



SUBMILLIMETER OBSERVATIONS OF MILLIMETER BRIGHT GALAXIES DISCOVERED BY THE SOUTH 

POLE TELESCOPE 

T. R. Greve, 1 J. D. Vieira, 2 A. Weiss, 3 J. E. Aguirre, 4 K. A. Aird, 5 M. L. N. Ashby, 6 B. A. Benson, 7 ' 8 
L. E. Bleem, 7 ' 9 C. M. Bradford, 10 M. Brodwin, 6 J. E. Carlstrom, 7 ' 8 ' 9 ' 11 ' 12 C L. Chang, 7 ' 8 ' 12 S. C. Chapman, 13 
T. M. Crawford, 711 C. de Breuck 14 T. de Haan, 15 M. A. Dobbs, 15 T. Downes, 2 C. D. Fassnacht, 16 G. Fazio, 6 
E. M. George, 17 M. Gladders, 711 A. H. Gonzalez, 18 N. W. Halverson, 19 Y. Hezaveh, 15 F. W. High, 711 
G. P. Holder, 15 W. L. Holzapfel, 17 S. Hoover, 7,8 J. D. Hrubes, 5 M. Johnson, 20 R. Keisler, 7 ' 9 L. Knox, 16 
A. T. Lee, 17 ' 21 E. M. Leitch, 7 ' 11 M. Lueker, 2 D. Luong-Van, 5 M. Malkan, 20 D. P. Marrone, 22 V. McIntyre, 23 
J. J. McMahon, 24 J. Mehl, 7 K. M. Menten, 3 S. S. Meyer, 7 ' 8 ' 9 ' 11 T. Montroy, 25 E. J. Murphy, 26 T. Natoli, 7 ' 8 
S. Padin, 2 ' 711 T. Plagge, 711 C. Pryke, 7,8,11 C. L. Reichardt, 17 A. Rest, 27 ' 28 M. Rosenman, 4 J. Ruel, 27 J. E. Ruhl, 25 
K. K. Schaffer, 7,29 K. Sharon, 7 L. Shaw, 30 E. Shirokoff, 2 ' 17 B. Stalder, 6,27 S. A. Stanford, 16 Z. Staniszewski, 2,25 
A. A. Stark, 6 K. Story, 7 ' 8 K. Vanderlinde, 15 W. Walsh, 6 N. Welikala, 31 R. Williamson 7,11 

Draft version March 8, 2013 

ABSTRACT 

We present APEX SABOCA 350 /jm and LABOCA 870 /im observations of 11 representative ex- 
amples of the rare, extremely bright (Si.4, nm > 15mJy), dust-dominated millimeter-selected galaxies 
recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with 
LABOCA with 40 < S^o^m < 130 mJy, approximately an order of magnitude higher than the canon- 
ical submillimeter galaxy (SMC) population. Six of the sources are also detected by SABOCA at 
> 3er, with the detections or upper limits providing a key constraint on the shape of the spectral en- 
ergy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified 
blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. 
These calibration samples inform the distribution of dust temperature for similar SMC populations, 
and this dust temperature prior allows us to derive photometric redshift estimates and far infrared 
luminosities for the sources. We find a median redshift of z — 3.0, higher than the z — 2.2 inferred for 
the normal SMG population. We also derive the apparent size of the sources from the temperature 
and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which 
supports the idea that these sources are gravitationally magnified by massive structures along the line 
of sight. 

Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: evolution — galaxies: 
starbursts — submillimeter — gravitational lensing 
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1. INTRODUCTION 
The first extragalactic surveys at submillimctcr 



SCUBA (Holland et al. 


1999) over a decade ago dSmail 


et al.|l997[|Hughes et al. 


1998; Barger et al. 1998 


) discov- 



ered a population of optf 
ies. Multi-wavelength follow-up campaigns subsequently 
established that the bright (Ssso^m ^3—15 mJy) sub- 
millimeter galaxies (SMGs) represent a significant pop- 
ulation of distant (z ~ 1 — 3) dust-enshrouded, far- 
infrared (FIR) luminous (Lfir ~ 10 13 L Q ) galaxies, 
in which large gas reservoirs (M(H2) ~ 10 11 M Q ) are 
being turned into st ars at a prodigious rate (SFR 
1000 M fi yr" ' 
20051. 



e.g. 



Neri et al. 2003 Chapman et al 



Ground-based extragalactic surveys conducted at 
wavelengths of 850 — 1200 /im have identified a few hun- 



dred sources in total sky area of ~ 4 deg (e.g., Cop 



pin et aL] 20061 Perera et al.||2008"[ |Austermann et aT 



2009| |Weitt et al.||2009J) The brightest sources found 
in these surveys have o'sso/im ^5 20mJy, and a source 
density of ~ ldeg~ 2 . A handful of brighter sources 
have been discovered in small area surveys by target- 
ing massive galaxy clusters and thereby taking advan- 
tage of gravitational magnification. For example, ultra- 
bright SMGs have been found behind the Bull et clus 
ter (IES 0657-56, S 87 n„.m = 48mJy, z = 2.79 ; IWilson 



et al [ 2008; Go nzalez et al.|l2010| ; [Johansson et"al.||201U| | 



and MACS J2135-010217 (SMM J 2135-0102, herTjifter 



called the "Eyelash", S 87n „. m = 106 mJy, z = 2.33; 


Swin- 


bank et al. 2010 

'l 1 1 — ' 1 ' 


[Ivison et al.|2010 


Danielson et al. 


2011 


)• 



found serendipitously (5*87-0 



Lestrade et al. 2010 



lOOmJy, z = 3 - 4; 
Ikarashi et al. 20111. These ob- 



um 



jects are important as gravitationally magnified windows 
on the obscured star formation process in the early uni- 
verse, but large samples of these rare objects cannot be 
obtained in surveys of a few square degrees. 

The landscape of submm surveys has changed dramati- 
cally in the last few years with the advent of observatories 
capable of mapping large areas of sky simultaneously in 



millimeter Telescope (BLAST - 
mapped oye r lOdeg 2 at 250,350 



Pascale et al. 


2008)) 


and 500 /Ltm 


(Devlin 



etaL] [2009 1 

vidmg photometric redshift estimates. The 10 m South 
Pole Telescope (SPT) surveyed 200 deg 2 to mJy depth 
at 1.4 and 2.0 mm in 2008, discovering a population of 
rare (~ 0.1dcg~ 2 ) and extremely bright (> 20mJy at 



1.4 mm) dusty galaxies (Vieira et al. 2010). Early re 



suits from the Science Demonstration Phase (SDP) ob- 
servations of the Herschel A strophysical Terahe rtz Large 



above, as "ultra-bright" SMGs. These sources appear 
more luminous than galaxies identified in smaller surveys 
(Ssso/jm — 3 — 15mJy); we refer to the fainter, un-lensed, 
population as "normal" SMGs. The apparent luminosity 
of the ultra-bright SMGs, and their excess relative to the 
expected high-flux number counts of normal SMGs, is a 
strong indication that a significant fraction are st rongly 



lensed with magnification factors of order 10 to 50 (Blain 
[19961 IBlain et al.|[l999l |Negrello et aLl|2007[ [Hezaveh k 
Hold er|201ip . 

Measuring this lensing amplification - using high- 
resolution follow-up imaging and detailed lens modeling 
- is one of the keys to determining the intrinsic properties 
of these ultra-bright SMGs and their place in the overall 
scheme of galaxy evolution. The other key, which we fo- 
cus on in this work, is determining their redshifts. Recent 
"blind" CO line searches, i.e., without prior optical/near 



IR sp ectroscopy, with the Zpectrometer ( Frayer et al. 
201 1\ on the Gr een Bank Telescope and Z-Spec (Brad- 
ford et al. 2009) on the Caltech Submillimeter Observa- 
tory yielded redshifts in the ra nge z ~ 1.6 — 3.0 for five 



ultra-bright H-ATLAS sources (|Frayer et al.||201l" Lupu 



et al.||2010[ ). At present, robust spectroscopic redshifts 
have been published for ten such ultra-bright SMGs - 
discovered either from the ground or with Herschel - al- 
though several more sources from H-ATLAS an d Her- 
MES have n ow been spectroscopically confirmed (Harris 
et al. 2012|. The bulk of these lie within the the ob- 
served reds hift range of normal, r adio-identified SMGs 



(z ~ 1 — 3, Chapman et al.||2005[ ), which may in part 
be due to the frequency coverage of Zpectrometer and 
Z-Spec favoring z < 4 CO detections. The redshift of 
the single z > 4 ultra-bright SMC published t o date was 



made using the IRAM PdBI WideX correlator ( Cox et al. 
2011[ ). A large overlap between the SPT and Herschel 



Area Survey (H-ATLAS - |Eales et al.||2010| and the 
Herschel Multi-tier ed Extragalactic Survey (HerMES - 

each surveying tens of deg 2 at ( |Spergel et al.||2003 l 



Oliver et al. 2010) 



sources is expected, but the longer SPT selection wave- 
length (1.4 mm) does predict a broader redshift range 
than the 350 — 500 /um-selected Herschel sources. The 
SPT may also be sensitive to a population of cooler 
sources that are invisible to the shorter- wavelength Her- 
schel selection. The large SPT survey area (2500 deg 2 , 
compared to the 600 and 380 deg 2 that H-ATLAS and 
HerMES cover) ensures the identification of the rarest 
and most highly magnified objects in the sky. 

In this paper we present the first 350 and 870 /im maps 
of a subset of 11 SPT sources. In p]the SABOCA and 
LA BOC A o bserv ations and data reduction are outlined. 
In ^3.1 and 3.2 we describe the sub m m m aps, source 
morphology, and fluxes, while in §3.3| and |3.4| we fit spec- 
tral energy distributions (SEDs), derive photometric red- 
shifts, FIR luminosities, and dust temperatures based on 
the SEDs. Finally, in p]we discuss the derived properties 
of SPT sources, including their redshifts, and the impli- 
cations these findings have. Throughout, we adopt a flat 
cosmology with = 0.27, = 0.73, and h = 0.71 



250, 350, and 500 /*m - al so identified analogous popu 



lations of b right sources (Negrello et al. 
et al.||201ip . 



2010 Conley 



in this paper, we denote SMGs with Ssso^m > 30 mJy, 
primarily discovered in the large area surveys described 

31 Insitut d'Astrophysique Spatialo, Batiment 121, Universite 
Paris-Sud XI & CNRS, 91405 Orsay Cedex, Prance 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. South Pole Telescope Selection 



The SPT ( |Carlstrom et al.||2011[ ), a 10-meter off-axis 
Gregorian design with a 1 deg^ field of view, has been 
surveying the mm-wave sky with unprecedented sensi- 
tivity and angular resolution since its commissioning in 
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2007. The SPT is located within 1 km of the geographical 
South Pole. At an altitude of 2800 m above sea level, the 
South Pole is one of the premier locations for mm-wave 
astronomy. The high altitude and low temperatures en- 
sure an atmosphere with low water-vapor content and 
excellent transparency. Meanwhile, the location at the 
Earth's rotational axis allows 24-hour access to the tar- 
get fields. 

The first receiver mounted on the SPT is a camera 
consisting of 840 transition-edge-sensor bolometers, op- 
timized for fine-scale anisotropy studies of the cosmic mi- 
crowave background (CMB) and the discovery of distant 
massive galaxy cluste rs through the thermal Sun yaev- 
Zel'dovich (SZ) effect |Sunyaev fc Zeldovich|[T972| .' The 
840 bolometers are split into six wedges each containing 
140 detectors. In 2008 (the season during which the SPT 
data in this work was taken), the array consisted of a sin- 
gle 3.2 mm wedge, three 2.0 mm wedges, and two 1.4 mm 
wedges. The 3.2 mm wedge did not produce science- 
quality data, but the 1.4 and 2.0 mm wedges performed 
to specification, resulting in r.m.s. survey depths of ap- 
proximately 1.3 mjy at 2.0 mm and 3.4 mJy at 1.4 mm. 
The 10-meter primary mirror of the SPT results in beam 
sizes (FWHM) of approximately 1.0' at 1.4 mm and 1.1' 
at 2.0 mm. The camera was upgraded to provide sensi- 
tivity in all three bands for observations starting in 2009. 
The 2500 deg 2 SPT-SZ survey was completed in all three 
bands in November 2011. 

The data r eduction pipeline ap plied to the SPT data is 



described in Vieira et al. (2010 hereafter V10). Briefly, 



the SPT time-ordered data from every working detector 
are calibrated and bandpass filtered, one or more com- 
mon modes are removed from the data of all the detectors 
on a given wedge, and the data from all detectors in a 
given wavelength band are co-added into a map using 
inverse-noise weighting. 

Sources are identified in SPT maps as described in V10. 
Objects are identified by first con volving the single-band 



map with a matched filter (e.g., Haehnelt & Tegmark 
1996 1 that de- weights noise and astrophysical signals on 
large scales and maximizes sensitivity to point-like ob- 
jects, then searching the filtered map for high-significance 
peaks. Sources are then cross matched between the 1.4 
and 2.0 mm catalogs, and their fl uxes are deboosted ac 



cording to the method detailed in Crawford et al. (2010 1 . 
Sources detected in SPT maps are classified as dust- 
dominatcd or synchrotron-dominated based on the ratio 
of their flux in the 1.4 mm and 2.0 mm bands. Approxi- 
mating the spectral behavior of sources between 1.4 and 
2.0mm as a power law, 5(A) oc A~ Q , we estimate the 
spectral index a for every source and use a — 1.66 as the 
dividing line between dust- and synchrotron-dominated 
populations, thereby removing all flat spectrum radio 
quasars from the sample (see V10 for details). We ap- 
ply an additional selection on the sample by impos- 
ing a cut on sources found in the Infrare d Astronomy 



Satel lite Faint-Source Catalog (IRAS-FSC, |Moshir et al. 
1992|. This corresponds to cutting out sources with 
> 200 mJy, which should remove any source at 
z < 1 from this sample. Absolute calibration for both the 
1.4 and 2.0 mm bands is derived from the CMB a nd the 
calibration u ncertainty is & 10%, as described in |Vieira 
elallpoTO). 



For the study presented here, 11 sources from 200 deg 2 
of 2008 survey data were imaged at 870 and 350 /im. 
Their 1.4 mm fluxes range from 17 to 40mJy and ap- 
proximate a flux limited sample. The sources are listed 
with their full source names in Table [T] throughout the 
paper we refer to them by their truncated coordinates 
(e.g., SPT-SJ233227-5358.5 becomes SPT 2332-53). 

2.2. APEX submm continuum imaging 

The submm observations presented in this paper were 
carried out at 350 and 870 /im with the Submillimeter 
Apex BOlometer CAmera (SABOCA) and the Large 
Apex BOlometer CAmera (LABOCA) at the Atacama 
Pathfinder Experiment (APEX) p^T he sources observed 
and their coordinates are given in Table [T] 

LABOCA is a 295-element bolometer array (Siringo et 
al. 2009) with an 11.4' field-of-view and a measured angu- 
lar resolution of 19.7" (FWHM). The center frequency of 
LABOCA is 345 GHz (870 /mi) with a passband FWHM 
of ~ 60 GHz. The measured noise performance for these 
observations was 60mJys 1 / 2 . 

SABOCA is a 39-element transition edge sensor (TES) 
bolometer array (Siringo et al. 2010) with a 1.5' field- 
of-view and a measured angular resolution of 7.8" 
(FWHM). The center frequency of SABOCA is 860 GHz 
(350 //m) with a passband FWHM of - 120 GHz. The 
measured noise performance for these observations was 
ISOmJys 1 / 2 . 

SABOCA and LABOCA observations of seven of our 
sources were carried out during Max-Planck observing 
time in May 2010 (PI: Wcifi). A further four sources were 
observed in August 2010 under ESO programs 086. A- 
1002A and 087.A-0968A (PI: Greve). This brings the 
total number of sources presented in this paper to 11. 
All SABOCA observations were obtained in excellent 
weather conditions (precipitable water vapor, PWV, less 
than 0.7mm). The LABOCA 850 /mi observations of all 
targets were done in good weather conditions (PWV < 
1.5 mm). 

Typical integration times for LABOCA were 1 h on- 
source resulting in an r.m.s. level of ~ 6mJybeam _1 , 
adequate to detect any of the SPT sources at z < 9, as- 
suming a point-like source (with respect to the LABOCA 
beam) and an SED similar to other SMGs. 

At the time of our observations, we had no informa- 
tion on source redshifts, SED shapes, or resolved struc- 
ture, all of which can dramatically vary the expected 
350 /im flux density for a given 1.4 mm flux density. The 
SABOCA observations therefore targeted a fixed r.m.s. of 
30mJybeam _1 (^3h on-source). At this noise level we 
expected a 5-er detection for a typical redshift of z ~ 3, 
assuming an unresolved point source. For sources that 
were not detected at this noise level we increased the in- 
tegration time up to 5 h on-source, yielding a r.m.s. of 
~ 20mJybeam~ 1 , which should detect a typical source 
out to approximately z ~ 4. 

Mapping was perfo rmed using the raster-spiral mode 
(Siringo et al. 2009[ ) for both bolometer arrays. This 
mode yields map sizes slightly larger than the field of 
view of the arrays (w 1.5' and 12' for SABOCA and 

32 APEX is a collaboration between the Max-Planck Institute 
fur Radioastronomie, the European Southern Observatory, and the 
Onsala Space Observatory. 
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LABOCA respectively). The SABOCA maps are thus 
well matched to the SPT resolution (FWHM ~ 1') and 
positional uncertainty of the SPT source (r.m.s. ~ 10"). 

Calibration was achieved through observations of 
Mars, Uranus, and Neptune as well as secondary calibra- 
tors and was found to be accurate within 10% and 25% 
at 870 and 350 /im, respectively. The atmospheric atten- 
uation was determined via skydips every ~ 1 — 2 hours as 
well as from independent data from the APEX radiome- 
ter which measures the line of sight water vapor column 
every minute. Focus settings were determined typically 
once per night and checked during sunrise. Pointing was 
checked on nearby quasars and found to be stable within 
3" r.m.s. 

The data were reduced using the Bolomet er Array 
analysis so ftware (BoA) reduction package (Schuller 
et al. |2010[ ). The time-ordered data undergo fiat field- 
ing, calibration, opacity correction, correlated noise re- 
moval on the full array as well as on groups of bolometers 
related by the wiring and in the electronics, flagging of 
unsuitable data (bad bolometers and/or data taken out- 
side reasonable telescope scanning velocity and acceler- 
ation limits), as well as de-spiking. Each reduced scan 
was then gridded into a spatial intensity and weighting 
map. Weights are calculated based on the r.m.s. of each 
time series contributing to a certain grid point in the 
map. Individual maps were co-added with inverse vari- 
ance weighting. The resulting map was used in a sec- 
ond iteration of the reduction to flag those parts of the 
time streams with significant source signal. This guar- 
antees that the source fluxes are not affected by filtering 
and baseline subtraction. For sources that remained un- 
detected at 350 /im we flagged a region with 15" radius 
centered on the 870 /im LABOCA position in the 350 /im 
time streams. 

2.3. Spitzer-IRAC imaging 

Mid-infrared Spitzer/Infrared Array Camera (IRAC) 
imaging was obtained on 2009 August 02 as part of a 
larger Spitzer program (PID 60194; PI Vieira) to follow 
up bright submm galaxies identified in the SPT survey. 
The on-target observations consisted of 36 x 100 s and 
12 x 30 s dithered exposures at 3.6 and 4.5 /im, respec- 
tively. A large dither pattern was used for the 3.6 fjxa 
exposures, and a medium dither pattern was used at 
4.5 /im. This scheme was designed to provide 3.6 /im 
imaging sufficiently deep to detect even very distant 
galaxies, while also providing a minimal level of sensi- 
tivity at 4.5 /im to furnish infrared colors for low-redshift 
sources nearby. 

The da ta were reduced f ollowing the methods de- 
scribed in Ashby et al.| ( |200"9| . The corrected Basic Cal- 
ibrated Data (cBCD) frames were modified individually 
to eliminate column pull-down artifacts, and treated to 
remove residual images arising from prior observations of 
bright sources. The resulting pre-processed frames were 
then mosaiced with standard outlier rejection tec hniques 
using MOPEX under the control of IRACproc (Schus- 
ter et al.] 2006 ) to create two mosaics covering the field. 
The final mosaics were generated with 0"86 pixels, i.e., 
pixels subtending a solid angle equal to half that of the 
native IRAC pixels. The field covered with at least ten 
3.6 /im exposures is roughly 6' wide, centered on the SPT 
position. In this work, we show only the 3.6 /im images. 



2.4. Spectroscopic Redshifts 

To date, spectroscopic redshifts have been obtained 
for two sources in this sample. SPT 0538— 50 has a 
redshift of z = 2.783, derived from ionized silicon 
(Si iv A1400A) emission in an optical spectrum with 
VLT/X-SHOOTER (284.A-5029; PLChapman) and car- 
bon monoxide CO(7-6) and CO(8-7) with APEX/Z- 
Spec (086.A-0793 and 087.A-0815; PI: De Breuck, 
and 086.F-9318 and 087.F-9320; PI: Greve). Multiple 
source images surrounding a low-redshift elliptical galaxy 
are resolved in high-resolution 890 /im imaging with the 
Submillimeter Array (SMA) (2009A-S076, PI: Marrone). 
The foreground lens has a measured spectroscopic red- 
shift of z = 0.443. A paper detailing these measurements 
and providing detailed characterization of this source (in- 
cluding lens modeling) is in preparation (Bothwell et al.). 

SPT 2332-53 has a redshift of z = 2.738, derived from 
Lya and ionized carbon (C iv A1549A) in an optical 
spectrum with VLT/FORS2 from program (285.A-5034; 
PI: Chapman) and carbon monoxide CO (7-6) with 
APEX/Z-Spec. The redshift of the foreground lens, 
which is a cluster of galaxies, is z = 0.403. A paper 
providing a full characterization of this source, including 
unambiguous evidence for it being lensed by the cluster 
of galaxies, is in preparation (Vieira et al.). 

3. RESULTS 

3.1. Submillimeter maps and source morphology 

Fig. nl shows Spitzer/IRAG 3.6 /im postage-stamp im- 
ages ofour 11 sources, with the SPT 1.4mm, LABOCA 
870 /xm and SABOCA 350 /im signal-to-noise ratio (S/N) 
contours overlaid. 

All sources are detected at 870 /i.m at S/N !> 6. For 
6 of the 11 sources, the S/N > 3 in the SABOCA 
350 /im images near the LABOCA positions (Fig. []}. 
A further source (SPT 2357— 51) is marginally detected 
at 350 /im as it shows a S/N ~ 3 emission peak at 
the LABOCA centroid position. Of the remaining four 
sources, two (SPT 2349-56 and SPT 2319-55) are tenta- 
tively detected (S/N ~ 2.5) in the SABOCA map, while 
two show no signs of a SABOCA detection. 

All sources, with the exception of SPT 0529— 54 and 
SPT 2332— 53, are consistent with being point sources in 
both the LABOCA and SABOCA maps. Two LABOCA 
sources are associated with SPT 2349— 56 (A and B, see 
Fig. [T]), each is consistent with a point source. Inserting 
point sources into the bolometer time-streams and sub- 
jecting them to the same reduction steps as detailed in 
|2.2| reproduces the observed maps very well. 

SPT 0529— 54 appears to be resolved at 350 /im, with 
S/N ^3 — 4 emission extending to the north and east. 
This three-component substructure is enclosed within 
the area of a LABOCA beam, and the source appears 
point-like at 870 /im. 

For SPT 2332-53 the LABOCA 870 /mi emission is 
clearly extended in the east-west direction, spanning 
~ V . At 350 /im the source resolves into three sepa- 
rate, unresolved blobs. We denote these sources A, B 
and C, going from east to west (see Fig. [lj. Two addi- 
tional sources (D and E) are seen in the LABOCA map 
to the north and south of the central east-west source. 
Subtracting point sources at the A, B, and C positions 
in the LABOCA map reveals an additional LABOCA 
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source, denoted F in Fig. [T] While D, E and F are sig- 
nificant (^ 5 — 6<t) at 870 /im, no associated emission is 
seen in the SABOCA map. For a detailed analysis of the 
properties of SPT 2332— 53 we refer to Vieira et at, in 
prep. 

3.2. Millimeter and submillimeter fluxes 

The LABOCA source coordinates, the SPT flux den- 
sities at 1.4 and 2.0 mm, and the LABOCA/SABOCA 
flux densities at 870 and 350 /im for our 11 sources are 
listed in Table [TJ The SPT flux densities have been de- 
termined using the method described in V10. We report 
the peak flux densities, as none of the sources is resolved 
in the SPT maps. Submm flux density determinations 
are discussed below. For all source fitting, the absolute 
calibration error was added in quadrature to the statis- 
tical error. 

We do not report the 2.0 mm flux of SPT 2332-53 since 



it coincides with a m assive cluster of galaxies (see Van 
derlinde et al. (2010)), and the SZ decrement of the clus 
ter at this wavelength overwhelms any emission (Vieira 
et al., in prep). In the case of SPT 2349— 56, we disre- 
gard the B component, and attribute the SPT 1.4 and 
2.0 mm fluxes to the component A, which coincides with 
the SPT centroid. 

The 870 /im LABOCA flux densities were determined 
by fitting a Gaussian with FWHM = 21.5", the resolu- 
tion of the smoothed maps, to the sources (except for 
SPT 2332— 53, see below). A constant baseline term was 
included in the fit in order to account for any background 
level. 

The 350 /im SABOCA fluxes were derived by sum- 
ming the flux within an aperture corresponding to the 
LABOCA beam solid angle (fi = 1.133 x FWHM 2 ), i.e. 

within a radius of r ap = y/Q/ir. A background level was 
estimated from a surrounding annular region (inner and 



outer radii of ri 



and r2 



,) and removed. 



The aperture was centered on the 870 /im LABOCA cen- 
troid. The flux uncertainties were estimated for each 
source by measuring the flux within the same aperture at 
the same position in 100 difference maps generated from 
individual observations, and measuring the variance of 
the resulting distribution. For the SABOCA detections, 
we compared fluxes obtained by centering the aperture 
on the SABOCA and LABOCA centroids, and found the 
differences to be < 10 %, well within the photometric un- 
certainties. 

For SPT 2332—53, which is resolved into three sources 
in the LABOCA map, we derived the total flux within 
an aperture of r = 44". This aperture encompasses com- 
ponents A, B, C, and F. The same aperture was used 
to measure the corresponding total flux at 350 /jm in the 
SABOCA map. The 1.4 mm SPT flux contains the com- 
bined contribution from all four components, and possi- 
bly also the D and E sources seen in the LABOCA maps, 
because of the large SPT beam. For subsequent photo- 
metric analysis we ignore any contribution from the po- 
sitions of D and E to the total flux because of their large 
distance from the detected source. The D and E sources 
are not consistent with being part of the same lensing 
configuration as A, B, C, and F. The D and E sources 
do not have obvious IRAC counterparts. We have inves- 
tigated the possibility that they are artifacts from the 



LABOCA data processing, but every test indicates that 
they are real. 

SPT 0529-54 is extended in the SABOCA map, so we 
determined its 350 /im flux using an aperture with radius 
25". 

3.3. Spectral energy distributions and redshifts 

The majority of our sources do not have spectroscopic 
redshifts, which prevents us from determining their dust 
temperatures and luminosities. Our data are not well- 
suited to photometric redshi ft measurements that rely 
on LIRG SE P templates (e.g. |Silva et al.|1998[ [Chary 
Elbaz|2001[ ), as we do not have data in the MIR and IN IK 
where these templates have strong spectral features. The 
SED of dusty galaxies is nearly featureless at the spec- 
tral resolution of bolometer cameras, thereby hampering 
simple photometric redshift techniques. The most promi- 
nent feature, the frequency of the SED peak, cannot pro- 
vide a redshift without prior information on the dust 
temperature, as these qu antities are degener ate through 
the ratio T d /(1 + z) (e.g. |Blain et al.| ( |2003| )). Nonethe- 
less, some success in predicting spectroscopic redshifts of 
SMGs have been achieved by maximum-likelihood tech- 
niques that employ libraries of SED templat es of local 
star bursts and ULIRGs in their analysis (e.g. Aretxaga 
etllp507 |. 

To overcome this degeneracy, we have developed a 
method that uses the distribution of SEDs for sources 
of known redshift taken from the literature. We assume 
a simple grey-body form for the SED to derive Td for 
sources in three samples: 1) unlensed sources at z > 1 
with 350 /xm imaging, 2) lensed sources at z > 1 with 
350 /im imaging, and 3) SPT sources with known red- 
shifts. 

For the unlensed sources, we draw our sample from: 



• 16 sources from Kovacs et al. ( 2006 1 , an 850 /im 
selected sample with 1.4 GHz radio counterparts, 
optical spectroscopic redshifts, and 350 /im imaging 
from CSO/SHARC-IlR 



• 11 sources from Kovacs et al. (2010), a z ~ 2 
Spitzer IRAC and MIPS selected sample with 
Spitzer/IRS MIR spectroscopy, IRAM/MAMBO 
1.2 mm imaging and 350 /im imaging from 
CSO/SHARC-II. 



• 18 sources from Chapman et al. (2010), a 1.4 GHz 
radio selected sample with 850 /im detections, op- 
tical spectroscopic redshifts, and 250, 350, and 
500 /im imaging from Herschel /SPIRE. 

13 optic ally and radio selected AGN from other 



studies (j Be nford et al | |1999| |Beelen et al.||2006, 
Wang et al. 2008 2010!) with optical spectroscopic 



redshifts, mm photometry, and 350 /im imaging 
from CSO/SHARC-II. 

For the lensed sources, we draw our sample from: 

• 5 well studied len sed sources from the literature, 
Blamlp99; 



(eg 



Bcnford et al. 1999 Barvai 



nis feTvr son 2002 ) including IRAS FT0T14+4724 



33 Some of these sources also have 1.1mm photome- 
try from CSO/BOLOCAM and/or 1.2 mm photometry from 
IRAM/MAMBO, which we use, when available. 
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Fig. 1.— Spitzer/IRAC 3.6 /Ltm images (grey scale; 1.8' X 1.8' for most, 2.2' X 2.2' and 2.6' X 2.6' in the cases of SPT 2332-53 and 
SPT 2349—56, respectively) of the 11 SPT sources presented here, with the corresponding SPT 1.4 mm S/N contours overlaid (red contours), 
shown at 3,5,7, etc. Also overlaid are the same S/N levels (including S/N = —3 shown as dashed contours) for the LABOCA 870 (im 
(green contours) and SABOCA 350 fim observations (blue contours). The regions shown are centered on the SPT centroid position. The 
LABOCA and SABOCA beam sizes are represented in the lower left and right corners of the top left panel, respectively. 
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APM 08279+5255, and H 1413+117 (the Clover 
leaf). 350 /jm ima ging for these sources comes from 
CSO/ SHARC-II ( |Benford et al.|1999[[BeeTen et al 



2006 ) and millimeter photometry from a variety of 
literature sources. 



• 2 strongly lcnscd sources discovered behind massive 
galaxy clusters (the Eyelash and the Bullet) and 
imag ed with SPIRE ( |Ivison et aL]|2010| |Rex et al. 



2010). 



7 lensed sources discovered serendipitousl y in large 
extragalactic -fferscfteZ/SPIRE surveys flNe grcllo 
et alj|2010| lOmont et al.||2011[ |Cox et al.||20TT 
Conley et al.||201ip . 



For the two SPT sources, we use the measured spectro- 
scopic redshifts described in §2.4| 

We fit each source in each sample with a blackbody 
law, modified with a spectral emissivity that varies phys- 
ically such that the dust opacity reaches unity at fre- 
quency v c (e.g., |Blain etaL|2003| ) : 

/„oc \\-^{vjv c f\B v {T A ) (1) 

Here, -B„(Ti) is the Planck function. We fix the spectral 
index of the emissivity to /? = 2.0, and cri tical frequency 
to v c ~ 3000 GHz (A c ~ 100 /mi) following [Draine| ( |2006| . 
From these fits we derive a distribution of dust temper- 
atures for each sample from the sum of the probabil- 
ity distributions of the individual sources (Fig. [2|. The 
SED fits exclude any data at wavelengths shorter than 
A bs < 250 /mi (rest wavelength ~ 50 /mi for the highest 
redshift sources) to ensure similarity between the cali- 
bration samples and the data in hand for SPT sources in 
this work. Moreover, such a constraint is appropriate to 
ensure the applicability of our simple model, which can- 
not adequately represent emission from hot components 
that may become apparent at shorter wavelengths. We 
find that our simple greybody fits all sources well over 
this wavelength range. We attribute the difference in the 
dust temperature distributions between the lensed and 
unlensed sources to the former being selected near the 
peak of their dust emission (e.g., by SPIRE at 350 /*m or, 
in the case of the bright AGN, by IRAS at 60/100 /mi), 
which will bias the selection towards hotter (and more 
luminous sources). The unlensed sources were selected 
from cither mm, radio or IRAC data, and we intention- 
ally marginalized over these three selection techniques to 
mitigate any bias. 

The distribution of Td, apparent Lfir, and z for the 
three samples can be seen in Fig. [3] The lensed sources 
have the highest apparent luminosity at a given Td, as 
expected from the gravitational magnification, which in- 
creases the solid angle they subtend. The lensed sources 
appear to be offset in the Tfir vs. Td space, indicative of 
them being randomly sampled from the background and 
made more luminous by gravitational magnification. The 
lensed sources pulled from the literature are significantly 
hotter than the unlensed population. This is presumably 
due to selection effects. 

The detection limits for the types of observations used 
in the three samples are also shown in the right panel 
of Fig. [3J for the SED model employed in this work and 
assuming a dust temperature of 35 K. The apparent lu- 



minosity threshold decreases at high redshift for SPT se- 
lection, owing to the steep rise in the submillimeter SED 
of these objects for (3 — 2. 

The observed submm/mm SEDs of the SPT sources 
are shown in Fig. [4] In all cases, the 350 /im data point 
falls below the power-law extrapolation of the longer- 
wavelength data, giving the photometric redshift tech- 
nique a spectral feature from which to impose its con- 
straint. To infer a photometric redshift for an SPT 
source, we compare our data to the modified blackbody 
SED model described above. We randomly draw a value 
of Td from one of our distributions and fit the redshift for 
the SED described by this temperature. Repeating this 
10 4 times, we generate a probability distribution for the 
source redshift that marginalizes over the Td prior, and 
adopt the median and standard deviation of this distri- 
bution as the estimate of the photometric redshift and 
its uncertainty. The unlensed sample has a median dust 
temperature of ~ 34 K, while the lensed sample has a 
median dust temperature of ~ 46 K. We adopt the distri- 
bution from the unlensed sources throughout this paper 
under the assumption that lensing is randomly sampling 
the underlying unlensed population of sources. This is 
a conservative choice, since of the three distributions, it 
results in the lowest redshifts because of the Ti/(1 + z) 
degeneracy. Also, the dust temperatures derived for the 
two sources with spectroscopic redshifts fall in the middle 
of the unlensed distribution. The photometric redshifts 
derived from this analysis are indicated in Fig. [4] and in 
Table El 

As a point of comparison, we also fit each SPT s ource 
with redshi fted SED template s of Arp220 M82 ( |Silva 
et al.||1998| and the Eyelash ( |Ivison et aL||2010[ ). The 
Eyelash SED is modeled using a two-component dust 
model comprising two modified blac kbodies with Ta = 3 
and 60 K and j3 = 2, as described in Ivison et al. ( 2010 ). 



We adopt the redshift for which a given template pro- 
vided the best fit to the data (global xl minimum) as 
the best estimate photometric redshift (for that tem- 
plate). The photometric redshifts and uncertainties ob- 
tained from these three SED templates are shown in 
Fig. |4j Without additional photometry shortward of 
350 /im it is difficult to discriminate between different 
SED templates such as Arp220 or M82. These galaxy 
template redshifts are consistent with the temperature- 
marginalized greybody values. 

For the two sources with spectroscopic redshifts 
(SPT 2332-53 and SPT 0538-50), we can also directly 
check the photometric redshift estimates. In both cases 
we recover the true redshift to within the statistical un- 
certainties. 

3.4. FIR luminosities, star formation rates, and dust 

masses 

Adopting the p hoto metric redshifts and SED fits de- 
scribed above in §3.3| we can derive far-infrared lumi- 
nosities and star formation rates, and dust masses. The 
fits are shown as black curves in Fig. [3J and are seen to 
provide a good match to the data. For the two sources 
with spectroscopic redshifts, we can break the Ti/(l + z) 
degeneracy, and derive the FIR luminosities, star forma- 
tion rates, and dust temperatures, which are listed in 
Table [2] For the remaining sources, we derive FIR lumi- 
nosities, star formation rates, and dust masses by fixing 
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Fig. 2. — A comparison of the dust temperature distributions 
derived from calibration samples of: (red) two spectroscopically 
confirmed SPT sources, (blue) all lensed high-z sources in the lit- 
erature with 350 (im detections and spectroscopic rcdshifts, and 
(green) unlensed 350 /im-detected SMGs with spectroscopic red- 
shifts. These distributions are used as priors on for our pho- 
tometric redshift technique, with the unlensed (green) being the 
primary choice for subsequent analyses. The error bars include the 
statistical and absolute calibration uncertainties. 

the dust temperature at the median of the distribution 
for the unlensed sources (T<j = 34K). 
We determine dust masses according to 



(1 + z)k Vi 



[B„ r (T d ) - B Ui {T C mb{z)T\ (2) 



where S Ua is the flux density at the observed frequency 
v Q = v T {\ + z)^ 1 (which we here set to 345 GHz, the 
rough central frequency of the LABOCA bandpass). Z?l 
is the luminosity distance, /i is the magnification fac- 
tor, and Tcmb(z) is the cosmic microwave background 
temperature at redshift z (which, in principle, has to 
be included, althoug h it only changes the dus t mass 

We 



at the - 2% level 
adopt K Vr /m 2 kg~ 1 



Papadopoulos et al. fl2000|)) 
0.045 x (t-y/250 G Hz^ 



( Hildebrand 

1983 Kruegel & Siebenmorgen 19941, where P = 2.U is 
the adopted dust emissivity index (see f 3.3 1 . The dust 
masses (before correcting downward for the unknown 
magnification) are in the range ~ 0.4—1.4 x 10 10 M 
(Table [2]). We note that significant uncertainties are as- 
sociateawith these estimates due to the uncertainties in 
the dust emissivity. 

The apparent FIR luminosity is determined from the 
modified b lack body SEDs, assuming Td = 34 K, as de- 
scribed in §3.3| The SED is integrated over the rest wave- 
length range 8 — 1000 /im. This yields apparent FIR lu- 
minosities in range Tfir ~ 1.8 — 6.5 x 10 13 Lq (Table[2|. 
To derive total s tar formation rates, w e use the following 
conversion from Murphy et al. (2011): 



SFR 



1.49 x 10~ 1( V 



Al, 



yr" 



ipm[8- 1000 /xm] 



(3) 



The resulting apparent star formation rates are in the 
range ~ 2.7 x 10 3 - 9.7 x 10 3 M Q yr" 1 (Table||. No cor- 
rection has been applied to Lfir or the star formation 



rate for the unknown gravitational magnification factor 
and we use the factor of /x to indicate that these val- 
ues are upper limits. Even if the lensing correction was 
known, we have assumed that AGN do not contribute 
significantly to the FIR luminosity, and the star forma- 
tion rates would therefore still be upper limits. 
Assuming a gas-to-dust mass ratio of ~ 100, consistent 



with what is found for normal SMGs (e.g., Santini et al. 
2010 1, we infer total apparent gas masses in the range 
0.4-1.4xl0 12 M o 



and gas-depletion time scales 
M gas /SFR ~ 100-220 Myr. This is consistent 



Mgas ^ 

of idcpi 

with the lower-limit gas-depletion time scales (idepl > 
40 Myr) inferred fro m CO observations of normal SMGs 
(Greve et al. 2005). Note that the gas-depletion time 
scale estimate is independent of the lensing amplification 
factor as it enters in both the estimate of total gas mass 
and the star formation rate. 

4. DISCUSSION 

4.1. FIR luminosities and dust temperatures 

In this section, we compare the derived properties of 
our SPT sources with those of other ultra-bright SMGs 
from the literature with 350 /im measurements and spec- 
troscopic redshifts of z > 1. This comparison is sum- 
marized in Fig. |3j which shows Tfir vs. z and Td vs. 
Lfir for the SPT sources and other ultra-bright SMGs, 
normal SMGs, and AGN. 

As can be seen in Fig. [31 where we plot Tfir vs. Td, 
the two SPT sources with spectroscopic redshifts (and 
thus well-determined dust temperatures) are significantly 
cooler than other lensed, ultra-bright sources with similar 
FIR luminosities. The dust temp eratures (determined 
using the SED fit described in { 3.3 1 of the strongly lensed 
sources from the literature are in the range Td = 34— 67 K 
with a median of 46 K, while the SPT dust tempera- 
tures have a mean Td = 35 K, see Table [2] and Fig. [2j 
The fact that the SPT sources appear to have lower dust 
temperatures than other ultra-bright SMGs with simi- 
lar FIR luminosities is likely due to the shorter wave- 
length (< 1 mm) selection of the latter, which will be 
biased towards warmer temperatures. Also, the ultra- 
bright SMGs (and AGN) from the literature were se- 
lected based on their extreme brightness, which probably 
implies a preference towards warmer dust temperatures. 
It is also possible that strong lensing preferentially selects 
towards compact and warmer objects. Finally, differen- 
tial magnification may h ave an effect on the measured 
source temperatures 



(e.g. 



Blain 1999), but it is beyond 



the scope of this paper to address this issue 

For SPT sources with only photometric redshifts, we 
infer apparent FIR luminosities (integrated from 8 — 
1000 /im) in the range - 2 - 6 x 10 13 L . In Fig. J| 
we plot the FIR luminosities as a function of redshift7 
along with those of the above mentioned lensed and un- 
lensed SMGs from the literature. The SPT sources have 
apparent FIR luminosities comparable to those of other 
lensed, ultra-bright SMGs, such as the SPIRE sources 
with spectroscopic redshifts, which have apparent FIR 
luminosities in the range 3.9— 7.8 x 10 13 Lq . In compar- 
ison, normal SMGs have an average intrinsic FIR lumi- 
n osity of 2.5 x 10 1 2 Lq (with a range ~ 10 n -3xl0 13 La-, 



Kovacs et al. ( 2006[ )). Thus, the apparent luminosi- 
ties of ultra- bright SMGs (including those discovered by 
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Fig. 3. — Left: Apparent Lpm vs. for all unlensed (green symbols) an d len sed (blue symbols) 350 /jm-detected SMGs and AGN 
(squares and diamonds, respectively) with spectroscopic redshifts at z > 1 (see j ]3.3| for references). Apparent (i.e. lensed) FIR luminosities 
vs. dust temperature for the two SPT sources with spectroscopic redshifts presented in this paper (red symbols). For all sources, the FIR 
luminosities and dust temperatures were derived by fitting a modified blackbody law with = 2.0 becoming optically thick at A c < 100 fj,m. 
Right: The FIR luminosity as a function of redshift for the same sources. The 9 SPT sources presented here which lack spectroscopic 
redshifts are also plotted (black dots with yellow error bars) with their ph otom etric redshift estimates and associated uncertainties. The 
uncertainty on Lfir given the assumed photometric redshift is described in j ]3.3| which assumes T^d = 34K - the median dust temperature 
of the unlensed population. The green, blue, and red dashed curves show the limiting FIR luminosity as a function of redshift corresponding 
to the effective SHARC-II 350 fim survey limit (r.m.s. ~ 5mjy), 3x the SPIRE 500 fim confusion limit (30mjy), and the 3 cr SPT 2.0mm 
survey limit (3.9 mjy), respectively, given the SED model of this paper and assuming a 35 K dust temperature. The SPT survey, with its 
longer wavelength selection, is more sensitive to sources at the highest redshifts (z > 5) than Herschel and the large survey area of SPT 
makes it sensitive to sources that are amongst the most rare and most highly magnified. 
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SPT) are typically an order of magnitude higher than 
those of the normal SMGs. To a large extent, this differ- 
ence is expected to reflect the lensing amplification of the 
ultra-bright SMGs. In the case of the Eyelash, where the 
amplification factor is well-known (/i = 32), the intrin- 
sic FIR luminosit y obtained after c orrecting for lensing 
is 2.3 x 10 12 Lq (Ivison et al. 2010), which is typical of 
normal SMGs. 

4.2. SPT sources are strongly lensed 

High-resolution follow-up observations of extremely 
bright SMGs in the literature have sh own the major- 
ity to be strongly len sed, dusty galaxies ( Swinbank et al. 
20101 iFuet al.|2012D. It seems reasonable to assume that 
FT - 



the SpT sources, given their similarity to the literature 
sources, are also highly lensed, although, we cannot rule 
out that some fraction of them are intrinsically extremely 
luminous. 

For some of the SPT sources, there is indeed un- 
equivocal evidence that they are strongly lensed. In 
the case of SPT 2332-53, the LABOCA source is re- 
solved into three distinct SABOCA sources - a morphol- 
ogy which strongly suggests we are dealing with a sin- 
gle source undergoing strong gravitational lensing. This 
is confirmed in optical/near-IR imaging of this source, 
where the SABOCA sources are clearly seen to line up 
with a giant arc (Vieira et al., in prep.). Similarly in 
SPT 0529-54, the LABOCA source is resolved into an 
extended SABOCA source with multiple peaks, again 
indicative of strong lensing. Finally, in at least one other 
case (SPT 0538— 50), we see structure in the submm 
emission consistent with gravitational lensing (Bothwell 
et al, in prep.). A lack of high resolution (sub)mm imag- 
ing prevents us from constructing accurate lens models 
and deriving the intrinsic properties of the remainder of 
the sources. The mass scale of the haloes responsible 
for lensing the SPT population (i.e. galaxy-scale lensing, 
cluster lensing, etc.) is still poorly determined, although 
theoretical modeling indicates that the majority of the 
source s should be lensed by ma ssive elliptical galaxies at 
z ~ 1 (H ezaveh &; Holder|2011 ). Future publications will 
address the lensing statistics tor t his s ample. 

The SED fitting performed in §3.3| provides a way to 
estimate the average magnification of the SPT sources. 
The luminosity and temperature of the sources are re- 
lated by a version of the Stefan-Boltzmann law, mod- 
ified to account for the frequency-dependent emissivity 
of the sources, so that Lfir = ^R 2 oT\, where we as- 
sume that the SPT sources emit as modified blackbodies 
isotropically from a surface of radius R. The effective 
Stefan-Boltzmann constant, a e g, is determined for our 
model by setting a e s /o~ equal to the ratio of the integral 
of the modified blackbody over a perfect Planck black- 
body. For our model, with the opacity reaching unity at 
A c = 100 fim and T d = 34 - 36 K (as for the two SPT 
SMGs with spectroscopic redshifts), a c ff/o~ varies from 
0.51 to 0.54. Using this relation, we solve for the appar- 
ent effective radii obtained for the SPT sources, and also 
for the lensed and unlensed SMGs and dusty AGN from 
the literature^] Assuming that the intrinsic sizes and 

34 For simplicity we have assumed that the modification of the 
Stefan-Boltzmann law is the same for all sources, given by the value 
for T(j = 34 K, the median for the unlensed population. It varies 



dust temperatures of the SPT sources are, on average, 
similar to those of the unlensed SMGs, we can ascribe 
the differences in the apparent radii to different average 
magnification factors. From Fig.[5]we find an average size 
ratio of (i?sPT)/(-RsMG,unicnscd) - 4.7, which translates 
into a magnification ratio of (hspt) ~ 22. We stress, 
that while the apparent size estimates will depend on 
the optically thick transition- wavelength, A c , all sources 
are fitted with the same SED model, and so changing A c 
will not alter the relative trend seen. 

For the SPT sources with photometric redshifts, the 
implied lensing magnification strongly depends on the 
assumed dust temperature. Sampling the distribution of 
temperatures associated with the unlensed catalog, we 
find a broad distribution of magnifications for the SPT 
photometric sources with 67% of the total samples in- 
dicating lensing magnification fi > 5. Fixing the dust 
temperature for all sources to be 33K < T<j < 35K, as 
was found for the two SPT sources with spectroscopic 
redshifts, we find the mean magnification of the SPT 
sources with photometric redshifts to be Jl = 11. In order 
to make the SPT sample of photometric redshift sources 
consistent with no lensing JL < 1, it is necessary for the 
mean dust temperature to be greater than > 80K. 

4.3. The redshifts of ultra-bright SMGs 

From the 11 SPT sources presented in this paper, we 
find a median redshift of z = 3.0 with a 68% range of 
z = 2.0 — 4.6 using the modified blackbody fits (see 



§3.3 ). If we instead use the Arp220, M82, and Eyelash 
SET3 templates, we find in all cases median redshifts of 
z = 3.1 — 3.2. We e mpha size that while the photometric 
redshifts derived in [3.3 should not be considered more 
reliable than within Az ~ 1 on a source-by-source basis, 
our analysis does indicate than when taken as a sam- 
ple (which includes two spectroscopic redshifts) the SPT 
sources do have a very high average redshift. The use 
of the unlensed sample of SMGs as the prior is con- 
servative in that it results in lower redshifts than would 
be found by using the lensed sample. The agreement 
between the modified black-body and template-derived 
redshift distributions, and the comparison to the two so 
SPT sources with spectroscopic redshifts suggests that 
this assumption is reasonable. 

The redshift distribution for the SPT sources is shown 
in Fig.|6]for all three choices of the calibration sample for 
the modified blackbody fitting technique. The redshifts 
for a sample of normal, radio-ide ntified, 850 Atm-selected 



et al. 2008 2011 Knudsen et al 



SMGs are also shown (z = 2.2, |Chapman et al.||2005 [) 
where we have updated the distribution to include five 
additional unlensed SMGs spectroscopically co nfirmed 
to reside in the r e dshift range z = 4.05 — 5.3 (Capak 

ICoppin et a 



201 



2010) as well as newly obtained reds 



2008 
lifts 



IBanerji et al 



of SMGs in the spectroscopic 'redshift desert 
0.8 — 1.5). If the SPT sources have similar dust tem- 
peratures as the unlensed SMG population at z > 1 with 
spectroscopic redshifts, then the FIR colors imply that 
the SPT sources are at systematically higher redshifts. 
For all choices of calibration sample (which provide the 

by a factor of ±50% over the 20 — 50 K range of T<j enclosing the 
bulk of the unlensed SMGs and asymptotes to unity at large T^, 
but factors of order unity are unimportant for this argument. 
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Fig. 4. — Observed submm/mm SEDs for the 11 sources presented here. The black symbols are SABOCA (350 fim), LABOCA (870 /jm) 
data, and SPT photometry (1.4 and 2.0 mm). The black curves are modified blackbody fits with fixed spectral index (/3 = 2.0 ). and 
Ta = 34K t aken from the median value of the unlensed sample of SMGs. Spectroscopic redshifts break the Td/(1 + z) degeneracy (Blain 
|et al.|2003[ l, and in those cases we quote the best-fit dust temperature. The red, blu e, and green curves represent the best-fit SEDs to the 
data based on the Arp220, M82, and Eyelash SED templates, respectively, (see j j3.3| , where the absolute scaling and the redshift has been 
allowed to vary (the latter is given in each panel). The greybody fit in black is the simplest curve which fits the data. Without additional 
photometry shortward of 350 /im it is difficult to discriminate between different SED templates such as Arp220 or M82. 
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Fig. 5. — The apparent radius (parametrized as <x {L-pm/T^) 1 ^ 2 
in units of kpc) of the two SPT sources with spectroscopic redshifts 
(red symbols). Other strongly lensed sources (SMGs and dusty 
AGN — blue symbols) from the literature are also shown, as are the 
normal, unlensed SMGs (green symbols). In the text we argue that 
the larger apparent radius of the SPT sources is due to gravitational 
lensing. 
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Fig. 6. — A comparison of our estimated redshift distributions 
of the SPT sources assuming a distribution derived from: (red) 
the two spectroscopically confirmed SPT sources, (blue) all lensed 
high-z sources in the literature with 350 fim detections and spec- 
troscopic redshifts, and (green) unlensed 350 /xm-detected SMGs 
with spectroscopic redshifts. The black solid histogram shows the 
normalized spectroscopic redshift distribution of normal u nlensed 
850 qm-selected S MGs bas ed on an updated v ersion of the |Chap^] 
|man et al.| | |2005| l sample ( |Banerji et al.|[2011| l , and including the 
hve spectroscopically confirmed normal z > 4 SMGs. Note that the 
black histogram are for individual spectroscopic redshifts, while the 
colored lines correspond to the probability distributions of photo- 
metric redshifts for the SPT sample. All curves have been normal- 
ized such that the integral is equal to unity. 

prior) , the SPT sources are a t higher median redshift 
p005|) 



than the |Chapman et al.| ( |2uua| ) sources. 

The SPT sources are hypothesized to be a sample of 
gravitationally magnified members of the normal SMG 
population, which suggests that they should be drawn 



from the same redshift distribution. The difference in 
their observed redshift distributions is likely due to a 
combination of selection effects. First, the |Chapman 



et al. ( |2005[ ) sample is strongly selected against high red- 
shift objects by virtue of the reliance on radio detec- 
tions to inform spectroscopic followup. Some correction 
is provided by the addition of newer sources found to 
be at higher redshift, but the selection remains impor- 
tant. For the SPT sample, we removed any source which 
is detected in the IRAS-FSC, thereby removing low red- 
shift unlensed interlopers and possibly a small number of 
lensed sources at z ~ 1. The longer- wavelength selection 
of the SPT sample also admits very high-redshift sources 
that would be harder to detect at 850 fxm as the SED 
peak approaches the bolometer passband. Finally, the 
strong lensing condition itself is expected to introduce a 
bias towards higher redshift objects. Theoretical work 
predicts that gravitational magnification by clusters and 
galaxies shifts the redshift distribution of an observed 
source population toward h igher redshifts than an equiv- 
alent field-selected sample jBroadhurst et al.|1995 Blain 



1996 Negrello et al.||2007 Jain & Limal 



5. SUMMARY 

We have presented SABOCA 350 urn and LABOCA 
870 /im maps and fluxes of 11 ultra-bright SMGs from 
the SPT survey, resulting in the first FIR/submm SEDs 
of SPT sources. Employing a novel method, which makes 
use of a compilation of SEDs of 58 sources with spec- 
troscopic redshifts at z > 1 and with 350 fxm photome- 
try, we have derived photometric redshift estimates for 
our sources. In the two cases where we have spectro- 
scopic redshifts, they agree well with the photometric 
redshifts from this analysis. The SPT redshift distri- 
bution, with a mean of z = 3.0 and an 68% confidence 
region z = 2.0—4.6, is found to be skewed towards signifi- 
cantly higher redshifts than the distribution observed for 
the normal un-lensed SMG population. We argue that 
the higher mean redshift of the SPT sources is due to 
the longer selection wavelength and the lensing selection- 
criterion. 

Using the median dust temperature measured from the 
unlensed population of sources from the literature, we 
fit modified blackbody laws to the observed SPT SEDs 
and estimated their apparent (i.e., lensing-uncorrected) 
FIR luminosities (£fir ~ 3.7 x 10 13 Lq), star forma- 
tion rates (SFR £ 2500 M yr _1 ), and dust masses 
(M d ~ 0.4 - 1.9 x 1O 1O M ). Based on the two SPT 
sources with spectroscopic redshifts, we derive dust tem- 
peratures (~ 35 K) similar to that of normal SMGs. 
From the same two sources we estimate a value of ~ 22 
for the gravitational magnification factor. Such a large 
magnification factor is in line with tentative lens models 
(Vieira et al. in p rep, Bothwell et al. in pr ep) and the- 
oretical modeling ( Hezaveh fc Holder] 2011 [ ). Once cor- 
rect ed for lensing, our fi ndings, as well as those of others 



(e.g. |Conley et alT][2011[ ) , suggest that the intrinsic prop- 
i-br 



erties of the ultra- bright SMGs will fall within the range 
of normal SMGs. Based on theoretical models and the 
measured projected area, it appears that the wide field 
coverage of the SPT finds the rarest, brightest objects, 
including those with the highest magnification relative to 
strongly lensed SMGs and AGN found in the literature. 
Due to the extraordinary boosting of the flux and 
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angular size caused by the gravitational lensing, ultra- 
bright SMGs are ideal targets for high-resolution follow- 
up observations with the Atacama Large Millimeter Ar- 
ray (ALMA), Herschel Space Observatory, and the Hub- 
ble Space Telescope. These observatories can access the 
dust properties, gas kinematics, and star formation con- 
ditions on sub-kpc scales in these systems. These ob- 
servations were based on sources selected from 200 deg 2 . 
SPT has now surveyed 2500 deg 2 and an extensive multi- 
wavelength campaign is on-going. 
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Submillimcter Observations of SPT Sources 



TABLE 1 
Flux Densities of SPT sources 



ID 




R.A. 






Dec. 








S350/jm 
mJy 


S870Min 

mJy 


Si. 4mm 

mJy 




S2.0mm 

mJy 


SPT-S 1051 259- 




05:12 


57 


g 


— 59 


35 


39 


(-; 


22"? -1- 95 


92 ± 12 


22.7 ± 


3 


9 


5.5 ±1.1 


SPT-S J052903- 


5436.6 


05:29 


03 





-54 


36 


33 


3 


409 ± 35 a 


112 ± 11 




A 
1 


8 


9.2 ± 1.3 


SPT-S J053250- 


5047.1 


05:32 


50 


9 


-50 


47 


10 





353 ± 47 


127 ± 10 


40.8 ± 


5 


2 


13.4 ± 1.4 


SPT-S J053816- 


5030.8 


05:38 


16 


5 


-50 


30 


52 


5 


336 ± 88 


125 ± 7 


29.7 ± 


4 


6 


8.5 ± 1.4 


SPT-S J055002- 


5356.6 


05:50 


00 


5 


-53 


56 


41 


2 


35 ± 40 b 


59 ± 10 


17.3 ± 


4 


3 


3.9 ± 1.0 


SPT-S J055138- 


5057.9 


05:51 


39 


2 


-50 


57 


59 


4 


165 ± 27 


76 ± 11 


26.7 ± 


4 


2 


5.0 ±0.9 


SPT-S J231921- 


5557.9 


23:19 


21 


5 


-55 


57 


57 


6 


50 ± 10 b 


41 ±5 


17.5 ± 


4 


2 


5.4 ± 1.2 


SPT-S J233227- 


5358.5 


23:32 


26 


5 


-53 


58 


39 


8 


425 ± 39 c 


150 ± ll c 


34.4 ± 


4 


7 


d 


-A 




23:32 


29 


7 


-53 


58 


38 


4 


148 ± 21 e 


51 ±7 










-B 




23:32 


27 


6 


-53 


58 


42 


8 


92 ± 21 c 


38 ±8 










-C 




23:32 


25 


8 


-53 


58 


38 





137 ± 21 


43 ±8 










-D 




23:32 


29 





-53 


57 


59 





65 ± 21 b 


22 ±4 










-E 




23:32 


29 


1 


-53 


59 


31 





39 ± 21 b 


21 ±5 










-F 




23:32 


27 


1 


-56 


58 


17 





17 ± 21 b 


18 ±4 










SPT-S J234942- 


5638.2 


23:49 


42 


3 


-56 


38 


41 


2 


94 ± 41 b 


82 ±8 


20.3 ± 


4 


5 


5.0 ± 1.2 


-A 




23:49 


43 


1 


-53 


58 


17 





5±41 b 


28 ±8 










SPT-S J235338- 


5010.2 


23:53 


39 


9 


-50 


10 


05 


2 


24 ± 21 b 


43 ±6 


19.9 ± 


4 


8 


5.8 ± 1.3 


SPT-S J235718- 


5153.7 


23:57 


17 


1 


-51 


53 


52 


1 


74 ± 39 


46 ± 7 


19.8 ± 


4 


6 


4.3 ±1.0 



Note. — SABOCA (350 /zm) and LABOCA (870 /jm) flux densities of the 11 SPT sources presented 
in this paper, along with their SPT (1.4 and 2.0 mm) flux densities. The positions are the 870 fim 
centroid positions, except for SPT 2332— 53-A, -B, and -C where the 350 /im centroid positions are 
used. Unless otherwise stated, all 350um fluxes were measured within an aperture corresponding to 
the LABOCA beam solid angle (see j ]3.2[ |. The errors given are statistical only. For the SED fitting 
analysis we add the absolute calibration uncertainty in quadrature to the statistical uncertainty. 
a The flux is integrated within a 25" aperture covering the extended emission. 

b No 350 fim emission peak at S/N > 3. The 3 50 fim flux was measured within an aperture corre- 
sponding to the LABOCA beam solid angle (see j ]3.2| l. 

c The flux is integrated within a 44" aperture covering the -A , -B , -C, and -F components. 
d No 2.0mm flux measurements available for this source, see §3.2| for details. 
e The flux is measured within a 13" aperture. 



TABLE 2 

R.EDSHIFTS AND MODELED FIR PROPERTIES OF SPT SOURCES 



Short ID a photo- z spec-z 6 L FIR [8 - 1000/mi] SFR T d M d 









xlO 13 L n~ x 


xlO 3 M yr _1 fi' 1 


[K] 


xl0 10 M Q fi 


SPT 0512-59 


2.5 ± 1.0 




3.4 ± 1.0 


5.1 ± 1.2 




1.0 ±0.5 


SPT 0529-54 


2.4 ± 1.1 




4.9 ± 1.1 


7.3 ±0.7 




1.3 ±0.3 


SPT 0532-50 


2.9 ± 1.2 




5.7 ± 1.2 


8.5 ± 1.6 




1.4 ±0.4 


SPT 0538-50 


2.5 ± 1.0 


2.783 


6.5 ± 1.6 


9.7 ± 1.5 


34 ± 3 


1.1 ± 0.2 


SPT 0550-53 


3.4 ± 1.4 




1.9 ±0.7 


2.8 ± 1.0 




0.6 ±0.5 


SPT 0551-50 


2.8 ± 1.2 




2.9 ± 0.6 


4.3 ±0.6 




0.8 ±0.2 


SPT 2319-55 


3.8 ± 1.6 




1.8 ± 0.5 


2.7 ±0.6 




0.4 ±0.1 


SPT 2332-53 


2.5± 1.1 


2.738 


6.5 ± 1.6 


9.7 ± 1.3 


36 ±3 


1.4 ±0.3 


SPT 2349-56 


3.1 ± 1.3 




2.7 ±0.7 


4.0 ± 1.0 




0.9 ±0.5 


SPT 2353-50 


4.5 ± 1.8 




1.9 ± 0.8 


2.8 ± 1.2 




0.4 ±0.2 


SPT 2357-51 


3.3 ± 1.4 




2.0 ± 0.7 


3.0 ±0.9 




0.5 ± 0.7 



Note. — Best estimates of dust temperatures, spectral indices and dust masses as derived from 
the modified blackbody fits (see j j3.3| . The FIR luminosities, star formation rates, and dust masses are 
apparent and have not been corrected for gravitational amplification, fi, and should therefore be considered 
strict upper limits. For the sources with only photometric redshifts, the errors of the FIR luminosities, 
star formation rates, and dust masses have been derived assuming T d = 34K 

a Shortened source names used throughout the text, with truncated coordinates. Sources are listed in the 
same order as Table [T] 

b Spectroscopic redshift derived from VLT spectroscopic observations in the restframe UV, and CO line 
observations with Z-Spec. For these sources we adopt the spectroscopic redshift when calculating the 
L Fm , SFR, T d , and M d . 



